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Abstract 

We introduce the notion of dissipative periodic lattice as an optical lattice with periodically 

distributed dissipative sites and argue that it allows to engineer unconventional Bose-Einstein su- 

C3 ■ perfluids with the complex-valued order parameter. We consider two examples, the one-dimensional 

dissipative optical lattice, where each third site is dissipative, and the dissipative honeycomb optical 
^\ 

\ lattice, where each dissipative lattice site neighbors three non-dissipated sites. The tight-binding 



approximation is employed, which allows one to obtain analytical results. In the one-dimensional 
case the condensate is driven to a coherent Bloch-like state with non-zero quasimomentum, which 
breaks the translational periodicity of the dissipative lattice. In the two-dimensional case the con- 
densate is driven to a zero quasimomentum Bloch-like state, which is a coherent superposition of 
four-site discrete vortices of alternating vorticity with the vortex centers located at the dissipative 
sites. 

PACS numbers: 03.75.Lm; 03.75.Nt 
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I. INTRODUCTION 



Bose-Einstein condensates (BEC) in the optical lattices have become an universal lab- 
oratory to study various quantum phenomena of the condensed matter physics [l-4|. A 
spectacular achievement in this direction is the demonstrated single-site addressability of 
the two-dimensional optical lattices js], 6|, i.e. the atoms can be controllably removed from 
the lattice sites selected at will. Another intriguing recent development is observation of the 
unconventional boson superfluidity in the higher bands of the optical lattice It-O], described 
by a complex-valued order parameter ^d, therefore, beyond the no-node theorem of R. 
Feynman (see for instance, the review [1Q[|). 

The purpose of the present paper is to demonstrate that the single-site addressability 
in the optical lattices allows one to generate the unconventional boson superfiuids. To this 
goal, one can use the controlled removal of atoms by the technique of Refs. js], ^ from a 
periodic sublattice of the optical lattice with BEC loaded in it. It is shown below that the 
condensate is driven to a coherent non-local Bloch-like state described by a complex-valued 
order parameter (in the two-dimensional case) or the order parameter with nodes (in the 
one-dimensional case). This is achieved by the emulation of a coherent non-local dissipation 
by a local one, similar as in the recently proposed scheme for generators of the non-classical 
states of photons 

The quantum state engineering based on the dynamics in the open quantum systems 
13 1, in particular in the optical lattices jl4|, is a well-known idea. As distinct from the 



previous proposal of the quantum state engineering in the optical lattices (see also the 



recent study 



15j), where the inter-site atomic currents are coupled to reservoirs, in our 



case the dissipation acts locally and incoherently on a periodic sublattice of the lattice sites, 
whereas the condensate is driven to a non-local coherent state. 

It is known that the action of dissipation in conjunction with the nonlinearity can increase 
coherence of the quantum state of BEC. For instance, the possibility to engineer the order 
parameter of BEC by a localized dissipative perturbation (i.e. by a dissipative defect) was 



shown in Ref . 



16| . It was also shown that the quantum coherence of a strongly interacting 



BEC loaded in the double-well trap subject to the phase noise and particle loss can be com- 
pletely restored by engineering the parameters of the system and controlling the dissipation 
rate Similar results were obtained with the optical lattices, where the particle loss at 
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the boundary acting together with the non 
and formation of the discrete breathers 



inearity resulted in restoration of the coherence 
18| . Whereas in the previous setups the action of 
a localized or boundary dissipation was considered, in the present proposal we consider the 
action of a periodic dissipation which profoundly affects the physics in the optical lattice. 
Therefore, to distinguish such a setup from the usual optical lattice, it will be called below 
the dissipative lattice. Indeed, a strong periodic dissipation effectively changes the lattice 
space group, it results in a bigger unit cell for the dissipative lattice as compared to the same 
lattice without the periodic dissipation. The dissipative lattice can have both the dissipative 
as well as the (almost) non-dissipative coherent modes, which are Fourier-like expansions 
over some superpositions of the local modes in each unit cell. The long-term state of BEC 
in the dissipative lattice is a lowest energy subset of the effective dark states, i.e. the states 
almost unaffected by the dissipation as the result of the quantum Zeno effect [lol, 
also the reviews 



21 



20j (see 



22|). Finally, in contrast to the previous works, a weak nonlinearity in 



our case plays only an auxiliary role (see the honeycomb lattice example below), which is to 
break the energy degeneracy of the effective dark states and select one particular state. 

The paper is arranged as follows. In section [IT] we consider the simple case of the one- 
dimensional periodic dissipative lattice, where the action of a periodically distributed dis- 
sipative sites is discussed. Then, in section llllt we consider a more involved case of a 
dissipative honeycomb optical lattice. Section HVl contains the discussion of the main results 
and a general perspective on the dissipative periodic lattices. 



II. THE ONE-DIMENSIONAL LATTICE WITH A PERIODIC SUBLATTICE OF 
DISSIPATIVE SITES 

Here we consider the simplest case of one-dimensional dissipative lattice to illustrate the 
main idea. Specifically, we concentrate on the one-dimensional lattice where each third 
lattice site is dissipative with the same dissipative rate, see Fig. [H The selected periodicity 
of the dissipative sites is the most dense one within the class of the dissipative sublattices 
which still allow for the non-dissipative coherent modes in the unit cell (see Fig. [U^a)). We 
consider the tight-binding approximation with the nearest neighbor tunneling and a weak 
nonlinearity as compared to the tunneling and dissipation rates. Then the standard boson 



3 



Hubbard Hamiltonian applies, which in the notations of Fig. [T]^a) reads 

H = -J ^{a|, o(a„ _ + an,+) + a^^+a^+i _ + h.c.} 

n 

+ ^ + (an,o)^an,0 + (l) 

n 

where the index n enumerates the unit cells of the dissipative lattice. Note that the dis- 
sipative lattice period is D = 3d, where d is the period of the original (non-dissipative) 
lattice. 

We assume that the atoms are removed with a constant and uniform rate F from the 
sublattice of sites as indicated in Fig. [H which can be realized, for example, by the technique 



of Refs. 



Then, the state of the system is given by the density matrix p satisfying the 



quantum master equation in the Lindblad form (see, for more details, Ref. 23| and Ref. 24 1 
for a general discussion) 

~[H,p]+Tj2^K,o]p, (2) 



dt h' 



n 



where the Lindblad term is defined as 'P[a]p = apa' — ^[a'^ap + pa'a]. 

The reason to use a larger unit cell of the lattice in Eq. ([T]) and ([2]) is that, for a 
strong dissipation rate, the population of the dissipative sites is locked to that of the non- 



dissipative ones (see Ref. 25| for further details) and can be adiabatically eliminated. In 
this way, the problem of finding the ground state of BEC in the dissipative lattice is reduced 
to that of a modified non-dissipative lattice, where the boson operators in each unit cell are 
some coherent modes over the non-dissipative sites (cn,± below). Our use of the introduced 
notion of the dissipative periodic lattice leads to a significant reduction of complexity of the 
analytical analysis. 

We focus on the limit of strong dissipation T ^ J /h and neglect the nonlinearity assuming 
it to be too weak to have contribution on the considered time scale. Specifically, we assume 
the condition J/{hr) ^ U{a^a)/J where (a^a) is the average atomic filling of the non- 



dissipative wells. By per: 



brming the adiabatic elimination of the dissipative sites (ctn.o) (foi' 



the details, consult Ref. 12| , section II) the master equation is reduced to that for the non- 
dissipative sites only (with the reduced density matrix p/j), which is in the same form as Eq. 
([2]) with, however, a reduced dissipation rate 7 and a reduced Hamiltonian if/j (the terms 
depending on the operators of the dissipative sites are thrown away). The coupling to the 
dissipative sites in Eq. ([1]) suggests to define the coherent modes in each unit cell as follows 
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FIG. 1: (Color online) (a) A schematic depiction of the one-dimensional lattice with a periodic 
sublattice of dissipative sites (open circles). The lattice unit cell selection is shown by vertical 
lines, where in each unit cell the lattice sites are enumerated by "— " , "0" and "+" (shown by the 
subscript to the local boson operator a). The consequential numbers on top enumerate the unit 
cells and the vector labelled by D at the bottom is the dissipative lattice period, (b) Schematic 
depiction of the Bloch mode (fL _. Here u = -^=. 

2 ' v2-L 



Cn,± = ('^n,+ =t a„ ._)/v^, then the tunneling part of the reduced Hamiltonian becomes 

Hr = J 5^{Cn,-Cn+l,- - ci^+Cn+l,+ + h.C.} 



4 

n 

J 



^{ci^+{Cn+l,- - Cn-1,-) + h.C.}. (3) 
n 

Since only mode c„,+ is coupled directly (with the coupling coefficient a/2 J) to the dissipative 
mode a„o in each unit cell, the reduced master equation reads (Ref. section II) 

n 

For the times staring from t ~ I/7 the state of the system is defined by the global coherent 
modes which are decoupled from the dissipative modes. Indeed, let us define the coherent 



5 



lattice modes (which can be called the dissipative Bloch modes) as follows 

^K± = ^2^e^ Cn,±, (5) 

where L is the total number of unit cells (i.e. L = N/3 where N is the number of lattice 
sites, see Fig. [1]). The unitary transformation ([5]) leaves the Lindblad term in the master 
equation (jl]) invariant, i.e. Eq. (jl]) has the same form in the Bloch basis 

^ = -l[Hn,pn\+lJ2'^[^n,^]PR. (6) 

n 

Eq. ([6]) defines a subspace of the dark states, which consists of the Bloch modes decoupled 
from the dissipative modes (pk,+- To find such modes we rewrite the Hamiltonian ([3]) in the 
Bloch basis 



+J ^ |« sin '^k,-fk,+ + h.c. I . 



(7) 



It is seen that the dark subspace contains two Bloch modes y^o,- and ifL _ (more precisely, 
the Bloch modes with S = {k — kd)/L <^ 1 , for kd = or k = L/2, have significantly reduced 
decay rate due to weak coupling to the dissipative modes, since the coupling coefficient reads 
J I sin [^k) I ^ J\6\). Moreover, since the dissipative lattice modes (fk,- have a negative 
kinetic energy, see Eq. ([7]), the lowest energy mode is the Bloch mode ifL _. Note that the 
long-term ground state qdl _ breaks the discrete translational symmetry of the dissipative 

2 ' 

lattice (the translation by D), since in terms of the local basis we have 



ifiL _ - 

V2L 

' n 



which is schematically depicted in Fig. [T]^b). 

Let us summarize the main conclusions of this simple example. The dissipative optical 
lattice with a strong dissipation rate is effectively equivalent to an effective non-dissipative 
lattice, which has unusual properties of the long-term ground state of BEC. In this particular 
case the effectively negative kinetic energy of the non-dissipative sublattice results in the 
Bloch state with a non-zero quasimomentum, which breaks the translational symmetry of 
the dissipative lattice. Below we consider a more intriguing example of two-dimensional 
lattice, where the sublattice of the dissipative sites induces a long-term state of BEC with 
a vortex-like phase distribution. 



III. THE HONEYCOMB OPTICAL LATTICE WITH A PERIODIC SUBLAT- 
TICE OF DISSIPATIVE SITES 



Let us now consider the lattice with the hexagonal symmetry, namely, the honeycomb 
optical lattice. The specific choice of the lattice has two reasons. First, in the tight-binding 
limit, the honeycomb lattice reduces to just four nearest-neighbor sites with equal tunneling 
amplitude, where three sites are arranged in the form of an equilateral triangle and one more 
cite is placed at its center, see Fig.[2|^a). From this perspective, by application of the external 
dissipation to the central cite in each such triangle one arranges for the three non-dissipated 
sites coupled to it, thus the next possible number as compared to the one-dimensional case 
considered in section HIl Second, within the class of two-dimensional lattices, the lattices with 
the hexagonal symmetry are known for their intriguing properties, such as the appearance of 
the Dirac cones (i.e. the Dirac Hamiltonian) in the Bloch band intersections and the related 
topological phase transitions |26^27J, including the phase transition from the conventional 
to unconventional superfluidity |9[. Thus, one naturally expects the dissipative honeycomb 
lattice to possess an interesting long-term ground state. 

The simplest honeycomb optical lattice, which was experimentally realized 



is created 



by intersection of three laser beams of of the same frequency and equal intensities. When 
the lasers are blue detuned from the atomic transition frequency, the lattice minima are 
at the vertices of honeycombs. In this case, the optical lattice potential reads ^(x) = 
Vo X]j=i '^os(bjx), where x = {x,y), bi 2 are the basis vectors of the reciprocal lattice, with 
1)2 being the vector bi rotated by vr/S, and bs = b2 — bi. We note that, interestingly, the 
following lattice 

V..8l/„=in^(5i^)=in^(if)s,„^(if). (9) 

shown in Fig. [2]^a), has also the honeycomb shape, but possesses much more pronounced 
minima of the lattice wells (it can be realized with the recent optical lattice technology 29| . 
see also Ref. 30|). The results obtained below, however, apply to any of the optical lattices 
reducible in the tight-binding limit to the honeycomb arrangement of the lattice wells, where 
in each equilateral triangle there are three non-dissipated sites located at the vertices and 
the dissipative site located at the center, i.e. as in Fig. [2l 

Taking into account the discussion of section [Tll we will use from the start the nomen- 
clature deflned by the dissipative lattice, see Fig. I^b). In particular, the lattice wells are 
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FIG. 2: (Color online) (a) Contour plot of the honeycomb optical lattice ([9]) (the lattice minima are 
shown) . The dissipative sites are located at the centers of the four-site elementary triangles (shown 
here by the division lines). The vectors Di^2 give the basis of the dissipative lattice periods, whereas 
Qi^2 connect the equivalent hexagons (shown by the thick lines in panel (b)). (b) Division of the 
dissipative lattice into the set of the equivalent hexagons used for representation of the Hamiltonian. 
The numbers at the vertices give the indices of the local boson operators corresponding to the non- 
dissipative wells, while the inner numbers enumerate the triangles inside each equivalent hexagon. 



enumerated as shown in Fig. |2l where each elementary triangle (shown by the thin lines in 
Fig. Et^a)) contains one dissipative lattice well at the center, to which we assign the boson 
annihilation operator ao, and three non-dissipative wells at the vertices, with the correspond- 
ing operators 01,2,3, as is schematically shown in Fig. [2]^b). It is seen that the dissipative 
lattice periods are given as Di = g^(2bi — b2) and D2 = 3^(2b2 — bi) (which are twice 
the original honeycomb lattice periods). It proves convenient to use the equivalent hexagons 
shown by the thick lines in Fig. EJ^b) as the unit cells of the dissipative lattice (see below). 
Therefore, we introduce the basis of the translations Qi = 2Di + D2 and Q2 = 2D2 + Di, 
Fig. [21^a), which translate the equivalent hexagons into each other without any change in the 
enumeration of the lattice wells. The annihilation operator at a given lattice site now is de- 
noted by aq,„,i, where / = 0, 1, 2, 3 is the lattice site index inside the triangle, n G {0, . . . , 5} 
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is the triangle index and q = giQi + q2Q2, Qj ^ {0, ±1, ±2, . . .}, is the hexagon index. 

Our division of the hexagonal lattice into the equivalent hexagons, as in Fig. [2t^b), and 
the specific enumeration scheme are dictated by the geometry of the strongly dissipative 
sublattice. Though the indices of the local operators assigned to the non-dissipative sites in 
each hexagon still can be distributed at will, the selected enumeration shown in Fig. Ulp) 
proves to be the most convenient. 

In the tight-binding limit with the nearest-neighbor tunneling confined to the lowest 
Bloch band, the boson Hubbard Hamiltonian reads 

(,) q,n,fc 

where the first sum is over the nearest neighbors and U is the nonlinear interaction strength 
proportional to the s-wave scattering length of BEG. 

We assume that the standard Markovian dissipation with the constant and spatially 
uniform rate F is applied to the lattice wells located at the centers of the elementary triangles 
of Fig. m Then the master equation reads 

% = ~[H.p\+TY,Y.V[a^,^,,]p. (11) 

q n=0 

In the limit of strong dissipation, i.e. 

the dissipative wells can be adiabatically eliminated, thus reducing the master equation f lTTj) 
to that for the non- d..pat,ve weUs only, S.ipp,,, the detaUs (see Ref. Q, seet.o. II)"* 

US write down the resulting reduced master equation for the density matrix pR describing 
the non-dissipative sites 

^ = + '^PR^ ^PR = ^ X] ^^[cq,„,3]PR, (13) 

q n=0 

where the reduced rate reads 7 = 12. P j (/i^F), H = Hr is the reduced Hamiltonian (the part 
of H from Eq. ([9]) dependent only on a^^n,! and aq^„ ^ with / = 1, 2, 3) and the local coherent 
basis used in the representation of Eq. f|T3|) is defined by the unitary transformation 

1 ^ 

Cq,„,/ = ^Y1 e'^'''«q,n,fc, 1=1,2, 3. (14) 

V-^ fc=l 



The reduced Hamiltonian Hr is obtained from the Hamiltonian (fTOj) by throwing away the 
terms with operators (aq,„,o and a^ci,n,o) belonging to the dissipative wells. Substituting from 
Eq. (I14p we obtain it in the form 



Hr = --jYl X^{^q,2.^^q,2s+l +^q,2.+l^*^q,2s+2 + h.C.} 
q ^s=0 

5 n 5 

+ E{^lnA^.+r„,„+3 + h.C.} (15) 



n=0 



q n=0 



where the vector notation ~(^l^^n = (Cq „ i, Cq_„ 25 <^li,n,3) is used. The coupling matrices are 
given as follows T^^^ = cos(|[Z — A;])e~'"'('"'^^/^, T* is its complex conjugate and Ai^k = 
|cos(^[/ — k]). The shift vector in the second sum of Eq. (fT5|) is defined as follows 
(see Fig. [2]) 

I"0 = Ql) I"l = Q2) r2 = Q2 — Ql, fn+S = —^n- 

The nonlinear interaction term in Eq. ( !T5|) is given in terms of the local coherent basis Cq,„^i, 
where the term H^^^ describes the nonlinear interaction in the elementary triangle (q, n). 
It has the following form (omitting the index (q, n) in Cq „j for simplicity) 

'-fc=l k<l 

+2(444 + 444 + 444 + h-c-)|- (le) 

We assume that the nonlinear interaction time scale is much larger than that of the 
reduced dissipation, i.e. tint ~ ^/{U{c^c)) ^ tdiss ~ 1/7; which implies the condition 

(note that the second condition in Eq. (fT2l) follows from condition (fT7|) and the first condition 
in Eq. f lT2|) ). Then the reduced Lindblad operator C of Eq. f|T3|) together with the reduced 
Hamiltonian (JT5l) define the long-time state |\E') of the condensate. First of all, |\E') must be 
the dark state of £: Cl"^) = 0. Indeed, this long-time state |\E') will be still subject to the 
weak dissipation coming from the nonlinear interaction in Eq. (|T6l) between the dissipative 
coherent modes Cq_„^3 and the coherent modes Cq^„ 1 and Cq^„_2 inside each elementary trian- 
gle. The form of the nonlinear dissipation can be found by a similar adiabatic elimination 
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procedure as above. However, since the rate of the reduced dissipation is proportional to 
the square of the interaction term (see Ref. jl2]), the rate of the nonhnear dissipation is 
proportional to the square of the interaction parameter, i.e. 7j„j ~ But then, under 

condition ( [T71) . the nonlinear dissipation can be neglected since its time scale is much larger 
even than the interaction time, we have l/7mi ~ tf^t/'^diss ^ tint, where tdiss is the reduced 
linear dissipation time scale and tint is the interaction time. 

Since the dissipative terms due to the nonlinear coupling can be neglected, the long- 
time state |\E') is a member of the dark states of C. The dark subspace is characterized by 
the set of lattice operators linearly decoupled from the dissipative operators Cq^n,3 entering 
the Lindblad operator C (|T3l) . To find these operators let us first diagonalize as much as 
possible of the quadratic part of Hamiltonian Hn ( |T5l) while keeping the Lindblad operator 
f|T3|) diagonal. This amounts to introducing the dissipative coherent lattice modes, i.e. the 
dissipative Bloch waves, by the following unitary transformation (note that there are exactly 
three types of different dissipative Bloch modes, which we combine in the vector notation) 

^ q n=0 

Here Af is the number of the equivalent hexagons in each of the two directions Qi^2 (see Fig. [2]) 
and p is the Bloch index, i.e. p = -^(piPi +P2P2) with pi 2 G {0,1,2, ... ,Af — 1}, while the 
reciprocal lattice periods Pi, 2 are defined by PjQj = 2TT6ij. They read Pi = |(bi — |b2) 
and P2 = |(b2 - ibi). 

The unitary transformation given by Eq. f|T8|) keeps the diagonal form of the Lindblad 
operator C in Eq. (fT3|) invariant, i.e. we have in the Bloch basis 

5 

^PR = l'^^1^[<fp,k,3]pR, (19) 
p k=0 

where the sum is now over the Bloch modes. Substituting Eq. ( ITSl) into the quadratic (i.e. 
tunneling) part of the Hamiltonian (fT5l) we obtain it in the form 



H 



(tunl) 
R 



t 



fci=0A;2=0 



(20) 
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where Sk,i is defined as 

^,,^(p) = e^('+-)l X^ef ('-)"cos (|[/ - m] + pr„) . (21) 

n=0 

One can easily check the invariance of the Hamiltonian Hji (ITSl) and the expression for its 
quadratic part ffj^*""') ( !20l) with respect to the dissipative lattice translation group, keeping 
in mind that a shift by Di 2 should be supplemented by a change in the operator indices. 
The latter amounts to the unitary transformation: Cq^n,k ~^ G'^^^Cci,n,k, where 6 = i.e. to 
the permutation of the enumeration of the wells inside each triangle according to the cycle 
(1,2,3) ^ (2,3,1), see Fig. ^b). 

Using the representation given by Eq. fl20|) one can find the Bloch indices (p, k) of the non- 
dissipated waves. Indeed, the operators (Pp,k,i and (Pp,k,2 which decouple from the dissipative 
operators v^p,fc,3, are solutions of the following equation (for / = 1, 2) 

^^Sk,M + Sk,,kM^i,3 = 0. (22) 

Since the first term in Eq. fl2^ is diagonal in (^1,^2), the matrix S'fc^^fc2(p) must be diagonal 
too. Expression ( 1211) gives only one case p = 0, i.e. the quasi-momentum must be zero. Then 
Eq. (!22|) reduces to a simple scalar equation with just two solutions for the second Bloch 
index: k = or k = 3. Therefore, the four Bloch waves which belong to the dark subspace of 
the linear dissipation correspond to the coherent operators v^o.o,; and <f 0.3,1 with / = 1, 2 in Eq. 
(HH]). They have the Bloch energies E{p = 0, A; = 0, /) = -2 J and E{p = 0,k = 3,l) = 2 J, 
respectively. Therefore the lowest energy is doubly degenerate. We note that such a dark 
subspace is an attribute of the dissipative periodic structure and would be impossible, for 
instance, with just one equivalent hexagon, since in this case the second term in Eq. f l22|) is 
absent. 

Due to the double degeneracy, the actual state of the system is defined, of course, by 
the nonlinear interaction part of the Hamiltonian ( IT5l) . In our approximation, given by 
Eqs. ( |T2l) and ( 1T7|) . we have U{a)a)/J <^ 1. Thus one can keep only the degenerate modes 
in the nonlinear term, which we denote as = v^co,! and = <^q,o,2- Projecting on the 
degenerate subspace, we obtain 

H^'""'^ = ^ {(^V)Vl + (^l)V^- + • (23) 
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FIG. 3: (Color online) Schematic depiction of the discrete vortex structure of the superfluid state 
|^+) (|24p . (a) The distribution of the vortices and anti- vortices in the elementary triangles of 
Fig. [21 (b) The distribution of the discrete vortex amplitude and phase between the wells (shown 
by contours) of the elementary triangles, with just two being shown, one containing the vortex and 
the other one - the anti- vortex. Here 9 = 



A simple analysis of Eq. ( I23l) shows that the interaction energy is minimized when all atoms 
are in either one of the two Bloch modes (p±. Each of the corresponding Bloch waves is a 
coherent superposition of the alternating discrete vortices and anti-vortices, each occupying 
just one of the elementary triangles, with the vortex centers being at the dissipative sites, 
as is easily seen from Eqs. (fT4|) and (fT8|) . see Fig. [31 In the explicit form 

5 

^ E E^^^^'l^'l.-.l) + ^^''\G^,n,2) + |Gq,.,3)}, (24) 
q n=0 

where 6 = 27t/3 and IGq.n,;) is the ground state in the respective lattice well. 

We note that though a macroscopic superposition of the two degenerate Bloch waves \'^±) 
is permitted by the model, one must recall that the discarded terms of the full Hamiltonian 
and the nonlinear dissipation will eventually destroy the macroscopic coherence, thus one 
of the two superfluid states will be spontaneously selected. For instance, only the coherent 
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superposition oc 1^1^+) — |^-) with only three bosons in the system is the only null state of 
the nonlinear dissipation. Thus, the long-time state of BEC in the dissipative honeycomb 
optical lattice of Fig. [2] has the complex- valued order parameter with the vortex-like phase 
distribution. 

IV. CONCLUSION 

In conclusion, we have introduced the notion of the dissipative optical lattice as the 
optical lattices with a sublattice of dissipative sites and have demonstrated its utility for 
engineering of the superfluid states of BEC with the complex-valued order parameter. The 
tight-binding approximation and the limit of a strong dissipation was used to allow for an 
analytical approach, since the numerical simulations would be difficult to carry out for the 
many-body open system on a lattice. We have found that the strong dissipation effectively 
changes the lattice space group, which results in a larger unit cell than that of the original 
optical lattice without the dissipation. Moreover, the base boson operators in the case of 
the dissipative lattice are linear superpositions of the local operators of the non-dissipative 
sites, what results in the ground state with a non-trivial complex-valued order parameter. 
For instance, in the dissipative honeycomb lattice the long-time ground state is a coherent 
superposition of the alternating discrete vortices and anti-vortices. 

We have considered just two examples, the simple one- dimensional case and the hon- 
eycomb lattice, where we have chosen some particular distributions of the dissipative sites 
which allow one to obtain the long-term ground state in explicit analytical form. However, 
similar results are expected to hold quite generally for the optical lattices with a sublattice 
of dissipative sites. Moreover, similar behavior of BEC is expected in dissipative optical lat- 
tices in the continuous limit (as opposed to the tight-biding limit), which will be considered 
in a future publication. 
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